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I’JTROD’u’CTIC!;  AND  NOTATION 


It  has  been  known  for  many  years  that  displayinr  the 
frequency  components  of  a  signal  as  a  function  of  tim.e  is 
a  valuable  signal  analysis  technique.  This  is  particularly 
true  for  nonstationary  signals,  such  as  speech,  v;hc5e 
frequency  structure  is  different  from,  moment  to  m.om.ent . 

In  the  case  of  speech,  the  spectrograph  has  been  widely 
used  to  plot  a  timie- frequency  representation  cf  a  spoken 
signal . 

In  this  paper  two  time-frequency  functions  will  be 
studied;  the  amibiguity  function  and  the  Ivigner  distribution 
These  functions  are  closely  related  and  yet  are  dissim.ilar 
enough  so  that  each  provides  a  valuable  perspective  to 
signal  analysis  problems.  Of  particular  interest  is  hov; 
these  functions  characterize  the  Heisenberg  uncertainty 
principle . 

§  1 .  Notation 

The  follovjing  notation,  definitions  and  theorem.s  will 
be  used  frequently. 

2 

a.  Given  f  €  L  CP),  the  Fourier  transform.  F  o! 
f,  sometimes  noted  as  F  =  f,  is  defined  as 

(1)  F(v)  =  fCtle"^"^^  dt . 

Here,  as  throughout  this  paper,  unless  othewise  noted 
the  limits  of  integration  are  assumed  to  be  from  -■»  to  +'» 

The  inverse  Fourier  transform  is  defined  as 


2 

f(t)  =  •“  J  F(v)  dv. 

^  2 

me  L  -norm  of  f  is 

'■f  =  (jl  f(t)!  ^ 

^2 

The  L  -norm  of  F  =  f  is 

"F  =  j|F(v)|^  dv)^^^ 

/\ 

Convolution  Theorem:  If  f  =  F  and  g  =  G  are 
Fourier  transform  pairs  and  h  =  fg  then  the 
Fourier  transform  of  h  is 

•)  f 

v)  =  1  F(r)G(v-r)dr  =  (F*G)(v). 

/n  j 

Plancherel's  Theorem:  Given  F  (  L^CIP.)  then 

2 

there  is  a  function  f  ^  L  GK)  such  that  F  = 
and  !!fii  =  '!Fl!. 

O 

The  Schwarz  inequality  for  L  -functions  is 
^  9  r  9  f  9 

ij  fgi  =  jifi  jur. 

^(t)  is  defined  to  be  the  Dirac  measure  having 
the  property  that 

Accession  Fn- 
NTIS  GtAA-r 


I  f(t)6(t)dt 


f(0)  . 


^♦i> 


IT.  AMBIGUITY  F’JUGTIO:: 

In  this  chapter  the  aiT'.blfuity  function  (AF)  will  be¬ 


st  uciea . 


of  its  more  useful  nronerties  will  be  stated 


and  rroved .  The  AF  v.’ill  be  shown  to  arise  naturally  frorr; 
the  output  cf  a  r.atche-l  filter  and  its  specific  application 
to  the  radar  probler.  will  be  discussed.  In  the  course  cf 
calculating  several  examples,  an  interesting  theorem  con¬ 
cerning  the  linear  transformation  of  the  coordinate  axes 
of  the  AF  will  be  stated  and  proved.  Finally,  the  AF 
will  be  related  to  signal  duration  and  shown  to  char.=  cterize 
the  Heisenberg  uncertainty  principle. 

The  AF  was  originally  introduced  in  l'^48  by  Ville 
in  the  context  of  a  general  signal  analysis  tool.  Applica¬ 
tion  of  the  AF  to  the  radar  problem  was  given  by  Woodward 
in  1950  [21].  Subsequently,  many  of  the  properties  and 
theorems  associated  with  the  AF  have  been  presented  by 
Siebert  [16]  and  Wilcox  [20]. 


§  2 .  Definitions 

The  AF  for  a  given  signal  f  has  several  comjnon 

definitions.  The  definition  m.ost  sim.ilar  to  the  Wigner 

2 

distribution  was  chosen  for  this  paper.  Given  f  €  L  OF) 
and  continuous ,  the  AF  X  of  f  is 


X(u,t) 


f(t  +^)f(t  -l)e 


T , _-lUt 


Although  the  definition  easily  generalizes  to  a  function  of 

2  ... 

two  continuous,  L  -functions,  such  a  definition  does  not 


iiave  a  sia.ple  rhyaical  interpretation  and  v.’ill  net  be:' 
studied.  If  f  is  not  continuous  it  v;ill  be  treated  as 
a  reneralieed  function. 

Other  definitions  are  sirdlar.  For  example,  in  the 
'cllcxinr  section,  the  AF  will  be  shown  to  be  the  output 
mciulation  of  a  matched  filter.  In  that  section  the 


ambi  r 

uity  function  wi 

11 

be  defined  as 

(2) 

0  (u 

,T)  = 

f 

J 

f(t)f(t-T)e 

X  is 

related 

to  0 

by 

the  formula. 

i 

(3) 

X  (U,T  ) 

= 

e  ^  0  ( u  ,T  ) 

In  radar  waveforrr.  design,  it  is  the  magnitude  of  the  AF 
v;hich  provides  the  measure  of  resolution  and  ambiguity  in 
f.  Thus,  in  this  context  definitions  (2.1)  and  (2.2)  are 
equivalent.  In  Papoulis’  text  the  AF  is  normalized  [12]. 
This  definition  is  not  selected  so  that  generalized  functions 
may  be  studied. 

Consider  now  the  expression  for  X  in  terms  of  the 
Fourier  transform  of  f.  First,  let  the  kernel  function 
Y  be  defined.  Given  f, 

(4)  Y(t,T)  =  f(t  +  p)f(t  -p)  . 

The  AF  of  f  is  the  Fourier  transform  of  y  for  fixed 
T .  That  is ,  Vt  f  H 


h.-; 


i 

I 


« 


h 


5 

(5) 

1 

X  (  U  ,T  )  =  j 

Y  ( t  ,T  )  e  dt 

■ 

This  e: 

\-nression  is  valid 

,  if  for  each  t 

€  P, 

Y  ) 

€  (IF.)  .  To  see  that  this  is  the 

case 

,  f  i  X  T  ^  T- 

and  use  Cauchy’-Schwarc  in  the  follo'winr; 

U)  (1 
J 

jY(t,T)!dt)"  = 

(j  1  f(t  1  ■ 

-y>id 

n 

t)  " 

< 

J 

f  of  0 

l|f|  l|fl  = 

1  J 

(f|f| 

J 

9 

Since 

f  €  L^CP,),  then 

y(*  ,T  )  e  L^(P.) 

for 

each  T  €  P 

and  (2.5)  is  valid.  Nov;  define  the  two-dimensional  Fourier 
transform.  T  of  the  function  y  ; 

(7)  r(u,v)  =  jj  Y(t,T)e  ^dtdT  . 

Calculate  this  integral, 

r(u,v)  =  II  f<t 

=  |(|  f(r)f(7;77e-""’^ 

If  f  =  F  then  f  (r-i )  F(-u).  Py  the  convolution 

theorem, 

,  .  f  ,  1  f  -pw - r  -isT  ^  -it(v-u/?), 

r(u,v)  =  (-^  F(-s)  e  F(u-s)ds)e  dr. 

Since  2iTS(v+a),  we  may  form.allv  write, 

r(u,v)  =  [  F(-5lF(u-s)5(s+v-u/2)ds . 


Hence 


r ( u , v) 


F(v+u/2)FTv-u72’5  . 


The  two-dimensional  Fourier  transform,  is  an  iterated 
integral,  consequently  (2.5)  and  (2.7)  may  be  combined  t 
rive , 


Y(t,T)  X(u,T)  r(u,v). 


By  Cauchy- Schwarz  we  see  that  for  fixed  u  €  IP. 


(1C) 


(|lr(u,v)|dv)^  5  (| 1 F(v) I ^dv)^ 


2  2 
Since  f  €  L  dP)  ,  by  Plancherel's  theorem,  F  €  L  dF.) 

and  we  conclude  that  for  every  u  €  IR  r(u,*)  €  L^dP). 

lie  may  then  write  from;  (2.9), 


(11)  X^(U,T) 


1  f  T/  UsivT, 

—  j  F(v  +-)r(v  --)e  dv, 


If  definition  (2.1)  is  used  to  define  the  AF  of  the 
Fourier  transform.  F,  we  have 


(12)  X_,(u,T) 


F(v  +i)F(v  -I)e‘^'^"'  dv. 


X^  is  related  to  X^  by  the  formula. 


X  ^  ( u ,  T  ) 


Thus  the  AF  of  the  function  f  is  a  scaled  and  rotate 


version  of  the  AF  for  its  Fourier  transform  F. 


AF 


in  Radar 


In  this  section  the  A.F  will  be  shewn  to  be  the 
complex  modulation  function  out  of  a  matclied  filter 
receiver.  This  will  help  interpret  the  tiiTie-frecuencv 
plane  on  v;hich  the  AF  is  defined  and  provide  a  foundation 
for  an  interpretation  of  the  function's  properties. 

A  i'adar  (an  acronym  for  R.Adio  Detection  And  Ranpinr) 
transmits  electromagnetic  energy  which  propagates  in  the 
atr.osphere .  Depending  on  the  reflective  properties  of  the 
objects  that  this  energy  contacts,  some  of  the  signal  is 
reflected  (or  reradiated)  back  to  the  receiving  system. 
Based  upon  the  difference  between  the  transmitted  waveform, 
and  the  reflected  v;aveform.,  the  radar  system  extracts 
information  about  the  target.  Of  interest  is  the  target's 
location  in  range  and  its  relative  motion  (radial  velocity) 

The  radar  vjill  be  modeled  as  a  matched  filter. 

Although  there  are  miay  ways  to  design  a  receiving  systemi, 
the  matched  filter  is  the  most  conunon  because  it  has  been 
shov7n  to  be  the  optimum  filter  for  a  very  large  class  of 
problem.s.  In  particular,  the  matched  filter  has  been 
shown  to  miaxim.ize  the  signal-to-noise  ratio  in  the  presence 
of  additive,  Gaussian  white  noise  [10].  Also,  it  has  been 
shown  to  maxim.ize  the  probability  of  detection  using  the 
m.axim.um  likeli?iood  criterion  [‘O.  Finally,  in  a  deter- 
m.inistic  sense,  given  a  known  signal  input,  the  m.atciied 
filter  has  been  shown  to  be  the  impulse  response  which 
miaxim.izes  the  output  at  a  riven  tim.e  with  respect  to  all 


ether  linear,  t irr:e- invariant  systems  [1?]. 

Trc  assumptions  are  made  in  this  deve  Icrr.ent  .  Fir 
tnc  sirnals  are  "narrowhand"  sirnals.  This  means  tnat 
river,  the  radar  siynal  is 

(1)  s(t)  =  a ( t )  cos  (at+C  ( t )  )  , 

then  a(t)  and  0(t)  vary  slowly  relative  to  the  hirh. 
carrier  frequency  w  .  (A  m;ore  precise  definition  cf 
narrov.'hand  may  be  found  in  reference  4).  This  perm.its 
the  cumbersome  expression  (3.1)  for  the  radar  signal  to 
written 

(2)  f (t)  =  c(t)e^^^, 

where  w  is  fixed,  the  complex  modulation  is 

(3)  c(t )  =  a(t)e^*^^^^  , 

and 

(4)  Ref(t)  =  s(t). 

(3.2)  is  the  complex  signal  representation  of  the  signa 
s(t)  and  is  discussed  further  in  Appendix  A. 

The  narrows  and  assum.pt  ion  is  also  necessary  to 
represent  the  doppler  effect  on  the  reflected  signal  as 
m.erely  a  shift  in  frequency.  The  doppler  aporoxim.at  ion 
develoned  in  ATcendix  .  Fortunately,  the  narrowband 
assuriotion  is  valid  for  almost  all  radar  applications, 
application  'where  such  an  assum.rtion  is  not  valid  is  in 


sonar  . 
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Secondly,  for  reasons  of  clarity,  the  transmitter  and 
the  target  from,  which  the  electromagnetic  energy  is 
reflected  are  riodeled  as  point  sources.  Modeling  the 
transmitter  as  a  point  source  is  accurate  since  the 
antenna  characteristics  do  not  effect  the  tine  delay  or 
doppler  shift  of  the  returned  signal.  Modeling  the  target 
as  a  point  source  ignores  relatively  small  effects  on  the 
returned  signal  which  do  not  serve  to  clarify  the  meaning 
of  the  ambiguity  function. 

To  begin,  define  the  transmitted  waveform 

(5)  f^(t)  =  c(t)e^‘^^. 


The  subscript  t  denotes  transmitted.  By  definition  of  a 
matched  filter,  the  im.pulse  response  of  a  linear,  time- 
invariant  filter  matched  to  the  transmitted  waveform  f^  is 

(6)  h(t)  =  f^(-t)  . 

If  f^  =  ,  then  the  transfer  function  of  the  matched 

filter  is 

H(u)  = 

The  radar  sends  out  f .  VJhat  returns  to  the  radar 
(assum.ing  a  target  is  present)  is  a  time  delayed,  frequency 
shifted  version  of  f .  The  time  delay  t  is  related  to 
the  range  of  the  target  R,  by 


t 


(7) 


T 


2R 

V 


V  is  the  velocity  of  propagation  of  the  signal  and  x  is 
the  time  it  takes  the  signal  to  make  the  trip  out  to  the 
target  and  back.  The  doppler  shift  u  is  related  to  the 
radial  velocity  v  of  the  target  by, 

(8)  u  = 

V 

(jj  is  the  fixed  carrier  frequency  in  radians  .  The  result¬ 
ing  reflected  signal,  when  it  enters  the  matched  filter,  is 

i(ai-u  )(t-T  ) 

(9)  ^r^^^  ”  c(t-T^)e  ^  ^  , 

Subscript  r  denotes  returned. 

The  radar  designer  wishes  to  maximize  the  probability 
of  detection,  hence,  would  like  to  match  the  receiver's 
filter  to  the  return  signal  f ^ .  But  the  parameters 
and  are  dependent  on  the  target  and  therefore  not 

known  a  priori .  Consequently,  the  designer  must  anticipate 
a  time  delay  x  and  a  freauency  shift  u  .  Then  the 
signal  to  which  the  filter  is  matched  is 

i(a)-u  )(t-T  ) 

(10)  f  (t)  =  c(t-T  )e  ^  ^  . 

m  m 

By  (3.6),  the  impulse  response  of  the  filter  matched  to 

f  is 
m 

-i(co-U_)  ( -t-T_) 

(11)  h  (t)  =  c(-t-T  )e  . 

m  m 

The  output  V  of  the  filter  h  with  the  input  f  is 

-  '  m  “  r' 


V  (  t  )  =  (  h  f  )  (  t  ) 

m  r 

r 

=  f  (x)h  (t-x)dx 
J  r  m 

r  i(u-U  )(x-T  )  -i(cj-u)(-t-T  +x) 

=  I  c(x-T  )c  ^  ^  C  (  -t-T  +x)  e 

J  r  TTi 

Let  n  =  X  -  T^, 

!•  i(u-u  )ri _  -i(oj-u  )(ri+T  -t-T  ) 

v(t)  =  |c(n)e  ^  c(n+T  -t-r  5” e  ^  ^  dn . 

J  r  m 

Let  T  =  t  +  T  -  T  and  u  =  u  -  u  ,  then 

m  r  r  m 

+  iT(u3-U  )  f 

(12)  y(t)  =  e  ^  j  • 

Compare  the  integral  on  the  r.h.-side  of  (3.12)  to  (2.2). 
We  see  that 

+  iT(to-u  ) 

(13)  y(t)  =  e  0(u,T )  . 


For  a  given  target,  and  t^  are  constants, 

0) ,  and  are  fixed  design  parameters  and  therefore 

the  AF ,  0(u,t),  is  a  function  of  time  (t  =  t+x  -x  ). 

Compare  (3.13)  to  (3.2)  and  we  see  that  y(t)  may  then  be 
interpreted  as  having  two  components;  the  carrier  and  the 
com.plex  modulation  function  0(u,x).  If  y(t)  is  viewed 
as  an  information  -  bearing  signal  we  might  expect  to 
extract  that  information  from  the  modulation  function 
0(u,x).  How  well  that  target  information  is  recovered  will 
depend  on  the  character  0(u,t). 
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Resolution  and  ar.biruity  are  two  measures  of  hov;  v;ell 
information  can  be  extracted  from  the  function  G.  Con¬ 
sider  first  resolution.  Given  design  parameters  u  and 
the  A"  at  a  specific  time  t  is  a  function  of  the 
tarret  parameters  u^  and  V.'ithout  loss  of  peneralitv 

we  can  study  the  0(U;T)  at  t  =  1.  Then  the  origin  of 
the  (u,T)-plane  represents  the  point  where  the  radar 

designer  has  precisely  anticipated  the  reflected  signal's 

2 

parameters.  Points  in  IR  away  from,  the  origin  represent 
mismiatches.  The  shape  of  the  AF  near  the  origin  will 
dictate  how  well  one  can  distinguish  the  actual  return 
from.  (0,0).  This  is  a  measure  of  the  resolution  of  f. 

Mext  consider  the  am.biguity  of  f.  It  will  be 
shewn  that  all  AFs  have  a  maximum  at  the  origin.  Should 
the  AF  have  peaks  away  from  the  origin,  then  for  a 
given  threshold  a,  the  set  E  =  {(u,t)  :  |0(u,t)|  >  a} 
may  be  disjoint  com.ponents  of  the  plane.  Ambiguity  then 
arises  in  determining  in  which  component  the  actual 
returned  param.eters  lie.  Both  concepts  will  be  discussed 
in  greater  detail  in  Section  6. 

§ 4  .  Properties 

1.  X  has  a  global  maximum  at  the  origin; 

(1)  |x(u,t)|  5  X(0,0),  V  (u,t)  CIF.^. 

Proof .  Let  g^(t)  =  f(t+^)e^'^^  and  =  f(t-  ^). 


Then  by  (1.6) 


(u,T)|'  =  I  g^(t)e2^t)dt|  5  |g^(t)l  dt  !f2(t)|^dt 

J  J  J 


:ub3tituti!'ir  for  and  we  obtain. 


(2)  |x(u,T  )!  ^  S  (f 


2.  ,2 

dt) 


(X (0,0)) 


The  result  nay  be  stated  as  |  X  ( u  ,t  )  |  t  '' f'  ,  since 


X  (0,0) 


°°  /, 

I f(t)I "  dt 

^  oo 


!I  2 


2  .  Synjnetry  . 


(  4  )  a  .  ) 


X (u,T )  =  X (-u,-T ) . 


Proof.  Note  that 


t)  =  f(t-^)f(t+~)  =  Y(t,T).  Then 


X(-U,-T 


-iut 


X(U,T  ) 


b.)  In  general  X  is  a  complex- valued  function.  If 
f  is  real  valued  and  either  f(t)  =  f(-t)  or 
f(-t)  =  -f(t)  then  X  is  real-valued.  That  is, 


X ( u  ,T )  =  X(u,T )  . 


Proof .  Only  the  case  of  odd  symmetric  f  will  be  verified 
The  case  of  even  symmetric  f  easily  follows.  Note  that 
f  real  and  odd  implies 

Y ( -t ,T )  =  f ( -t  +  ^) f ( -t  -  ^)  =  f(t+^)f(t-^)  =  Y(t,T). 

Then 


X  ( u ,  t") 


^  -  -  -  -  -  ^  -  —  1  n't" 

I  Y(t,T)e  dt  =  J  Y(-t,'i)e  dt  =  X(u,t) 


Translations  of  f(t)  and  F(u)). 


a.)  Let  T  be  the  translation  operator,  that  i 

3. 

T^f(t)  =  f(t-a). 

O. 

T  'n  0  r. 


(C)  X^(u,t). 

a 

Proof .  The  result  follows  from  the  fact  that 
'.T-  f  ( t  ,T  )  =  Y^(  t-a  ,T  )  . 


b.)  If  T^F(v)  =  F(v-b)  and  g(t)  •*-->  F(v-b)  th- 


(7) 


X^(U,T) 


ibT  ,  V 
e  X^(u,t) 


ibt 


Proof .  It  is  known  that  e  f(t)  F(v-b)  .  Hence 

-ib(  t-^)- 


Yg(t,T) 


ib(t  +-:r) 
e  f  (t  + 


ibT  . 

e  Y^(t,T) 


The  result  follows  by  the  definition  (2.1) 


4.  Modulations  of  f  and  F  by  simple  sinusoids. 
Let  M^  be  the  modulation  operator  defined  as 


M  f(t) 

=  f(t), 

a 

and 

:',F(v) 

=  r(v). 

(?) 


Then 


f(u,T) 


iai 

e 


X  ^( u , T  )  . 


Furtherrr.ore ,  if  r  .M,  F  then 

D 

(9)  X  (u,t)  =  X.(u,t) . 

Froof .  The  proof  of  (4.S)  is  the  same  as  that  for  (4.7), 
anf  the  proof  of  (4.9)  is  the  same  as  that  for  (4.6). 

Concise  statements  may  also  be  made  for  the  AFs  of 
Fourier  transforms.  For  example,  property  3b  m.ay  be 
written  as 

(10)  p(u,T)  =  e‘^^^  X^(u,T )  . 

Proof.  Let  g  *-■*  T  F,  then  from  (2.13), 

s. 

X„  r.(u,T  )  =  2ttx  (t  ,-u)  . 

1  i  g 

a  ^ 

From.  (4.7)  and  again  (2.13)  we  conclude, 

X„  p(u,T)  = 
a 

VJe  conclude  from  equations  (4.6)  and  (4.7)  that  the  modulus 
of  the  AF  is  invariant  to  translations  of  f  or  F. 

This  is  a  significant  feature  of  AFs. 

5.  Multiplication  of  functions. 

The  previous  results  can  be  generalized  to  modulation 


2TTe'^^'*  Xp(T  ,-u) 


-lau  ,  , 

e  Xp(u ,T ) . 


of  f  by  any  function  g 


(11)  X  .  ( U , T  ) 


2-  J  X^(r,T  )x^(u-r,T  )dr. 


Fr£o_f.  liotice  that  =  y  t  ,t  )y  t  ,t  )  then 


x^(u,t)  =  y^(t,T 

J 


)y  (t,T  dt . 

o 


Held  T  fixed  and  use  the  convolution  theorem  and  the 
result  follows . 


6.  Linear  filtering. 


h(t)  =  (f*g)(t), 


where  g  is  the  impulse  response  of  the  linear,  time- 
invariant  filter.  Then 


(12)  X^(u,t) 


X^(u,r)x^(u,T-r)dr 


Proof.  It  is  known  that  with  h(t)  so  defined, 


H(w)  =  F(w)G(a)), 


where  f  =  F  and  g  =  G.  Use  (2.11)  to  write, 


X^(u,T)  =  ^  j  H(v+H)  h(v-H)  d- 


-1-  T*/  _i_^\  T^/  ivx 

—  F(v+-)  F(v-~)G(v+-)  G(v-tr)  e 


The  result  follows  from  the  convolution  theorem,  and  ( 


and  (2.9)  which  state  that  for  fixed  u  f  IF 


>;  ^  ( u ,  T  ) 


Y  F(v  + -)  F(v  -  p 


an'.: 


X^(u,t)  y  G(v  -  y) 


7 .  Dilation . 

Let  g(t)  =  f(at),  for  5  oie  a  €  IR, 


(13)  then  X  ( u  ,t  ) 
F 


a !  X-  ( -• ,  aT  )  . 
^  a 


^roof .  ?y  definition. 


X  (u  ,T  ) 
F 


.c/  .4.  j.  X  jr  /  .  ai  X  -lut  , 
f(at  +-7r)  f  (at  — ^)e  dt 

/  2 


For  a  >  0  we  conclude  that 


-lU- 


X  (u,t)  =  I  f(r+^)  f(r-^)e  ^  -  dr 

£  J  2  2  a 


~  X.(-,  at) . 
a  fa 


For  a  <  0  we  conclude  that. 


X^(u,t ) 


So  Va  6  IR, 


X  (U,T) 
£ 


--T-  X  (-,  aa) 
a  fa 


8.  Invertibility  and  Uniqueness. 

Given  an  ambiguity  function  X,  the  generating 
function  f  may  be  uniquely  recovered  to  within  a  multi¬ 
plicative  constant  c  €  (C  such  that  jc]  =  1. 


.fciifc  ti»'i  »  ^  I 


Proof .  Assume  two  functions  f  and  g  generate  the  car. 
AF  so  that 

=  X  . 

Then  by  the  uniqueness  of  the  Fourier  transform  we  can 
conclude  that  =  Y  •  That  is 

j.  £ 

f(t+^)  f(t-^)  =  g(t+^)  g(t-^). 

z  ^  z  z 

t  +  “  and  t^  =  t  -  Then  (4.14)  is, 

g(t^)  "  ^ 

Hence,  \/t  €  IR  ,  g  =  cf  and 

f  _  cf 

p  -c  ”  _  * 

Cx  ^ 

1  -  1  I  I  2 

-  =  c  so  1  =  I  ci  . 

This  raises  the  question  as  to  how  one  identifies  an 
ambiguity  function.  That  is,  what  are  the  necessary  and 
sufficient  conditions  for  a  function  X(u,t)  to  be  an  A 
This  is  of  particular  interest  in  radar  waveform  design 
vjhere  one  would  like  to  establish  the  AF  so  that  X 
displays  the  desired  resolution  and  ambiguity  character¬ 
istics.  After  inverting  X.  the  designer  would  have  f, 
the  suitable  radar  signal.  Much  work  has  been  devoted  to 
this  effort  but  the  best  that  can  be  said  about  the 


(14) 

Let  t^  = 
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sufficient  conditions  for  X  to  be  an  anbifuity  function 
is  that  when  the  inverse  Fourier  transform,  is  -Derforiried  on 
X(u,t),  with  respect  to  the  u  variable,  the  result  is 
the  factored  form  of  the  generating  function  Y(t,T).  To 
be  TiOre  precise: 

9.  X(u,t)  is  an  amhiguity  function  if  and  only  if 

_  -I  (•  i  ^(t  +t  ) 

(15)  f(t^)  t  ( t  „ )  =  1  X(u,t2-t^)e  ^  ^  du . 


Proof .  The  necessary  and  sufficient  conditions  are  just 

a  restatement  of  the  definition.  Let  ''"i  ~  ^ 

t^  =  t  +  "j,  and  rely  on  the  uniqueness  of  Fourier  transforms. 

10.  The  squared-magnitude  of  X  has  the  unusual  property 
of  being,  after  a  coordinate  transform,  self-reciprocal 
in  the  two  dimensional  Fourier  transform.  That  is, 


(16) 


-i(uv+tT  ) 
e 


dud  T 


2m  I X(t  ,-v) 1 ^ . 


Proof .  Use  definition  (2.11)  for  x  term.s  of  F  =  f, 

and  property  (4.4)  to  obtain. 


(17)  x(u,t)  =  X(-u,-t) 


_L 

2m 


F(v  -p 


F(v  +^) 


-I 

e  dv 


( 


I 


i 


Use  this  expression  in  the  expansion  of  (4.16)  which  follows. 


^  -  M  -  ^-i(uv+T  t)  . 

X  (  ^  i,  /  I  0  U  li  u  i 


|j(  f(x  +])*(>:  -'^-)c  p  j  r(y-y)r(y+y) 

- i ( UV+  Tt )  .  . 

X  e  cue ■"  . 


u  V  - 1  ■/T 


Vee  the  chance  or  vai"ia^  les  r  =  >:  +  -  and 


f  f  !x(u,t  )! 

J 


2  -i ( uv+tT  ) 
e 


f(r)F(5)(  f  (r-T  )  e  F(s  +  u)e  ^ du )  drds 


rhe  first  inner  integral  is. 


I  -Tf - T  -iT(s  +  t) 

J  f(r-T)  e  < 


^-ir(s+t)  p(5+^)_ 


The  second  inner  integral  is 


1  f  w — T  ^  -iu(r+v)  , 
— -  F(s  +  u)e  du 

/"T 


eis(r*v)y;^rr7T 


Substituting  these  back  into  (4.18)  and  rearranging  terms 
we  see  that 


j  I X(u,T ) 1 


2  -Kuv+tx),  , 
e  dudx 


X^(t,  -v)  Xp(v.  t) 


The  result  follows  from  the  formula  (2.13). 


11.  The  radar  uncertainty  principle. 

2  2 

Integrating  |x|  over  all  of  IR  we  have  the 
interesting  result. 


X 


Proof .  This  IS  a  special  case  of  (4.16)  with  v  =  t  =  0  . 

This  property  can  be  interpreted  as  the  "conservation 
cf  ariii  raity  property."  It  says  that  the  best  a  radar 
de-'irner  can  do,  riven  a  specific  enerrv  constraint,  is  tc 
snifi  the  ar.biruities  inherent  in  f  to  the  unused  parts 
of  the  (u,T)-rlane.  Another  interpretation  is  that  the 
ar.ount  cf  ar.bipuity  in  f  is  invariant  over  the  class  of 

r\ 

functions  whose  b  -nor.t.  is  the  sane.  This  will  be  further 
discussed  in  Section  6. 


§  5  .  Examples  of _ . 

In  this  section  examples  of  AFs  will  be  calculated. 
Each  illustrates  some  of  the  properties  of  Section  4 .  The 
sirnificance  cf  some  examples  vjill  be  discussed  more  fully 
in  Sections  6  and  7. 


Example  1. 


Let  f(t)  =  ?j(t) 


1  for  t  <  T 


0  elsewhere. 


Then  by  (2.4),  y^TiT)  =  p.p(t+Y)Pj(t-Y).  Y(t,T)  is  unity 
inside  the  rhombus  below  and  zero  elsewhere. 


Firure  1 


7 

i 

" . 


\(u,t)  =  I 


( t  ,T  )e  dt  . 


J 


’  T  ^  [-2T,  2T]  tihe  lirrits  of  inte,rT''atiQn  arc 
+!t|/:  tc  t=T-  1x1/2.  Thus 


\  (  u  ,  T  ) 


i -T+l T 


_ut 
e  dt 


1x1/2 


ind  we  conclude 


1)  >(u,t) 


'  2  I  -  I 

■  -£ir.[u(T  -  for 


0 


.t|  -  2T 
for  1x1  ^  2T . 


lotice  that  for  x  =  0  and  u  =  0  V7e  have 


X(u,0) 


-sin  uT 
u 


for  all  u  €  IR. 


X  (  0  ,  t ) 


2T  - 
0 


T  1  for  I T 1  2  2T 
elsewhere . 


hese  are  sketched  in  Fipure  2 


Firure  2 


Exarr.ple  2.  Let  f(t)  =  e  p^(t). 

The  using  property  4  and  the  previous  example. 


(  2  )  X  ( u  ,t  ) 


-  Ein[u(T  - -L^)  ]  for  |t1  S  2T 
u  2  '  ' 


for  I  T  I  S  2" 


Example  3.  Let  f(t)  =  for  all  t  C  IF. 

Then,  formally. 


vc,,  ^  -  I  ^ia(t+T/2)  -ia(t-T/2)  -iut 

(3)X(u,T)  =  |e  e  e  dt 

J 


Recalling  1  2:15  (v),  (5.3)  is 


X  (  u  ,  t ) 


2TTe  6  (u)  ,  V  T  €  IF. 


This  means  that  the  AT  of  f  concentrates  all  of  its 
mass  on  the  u  =  0  axis.  Hence,  f  has  perfect  resolution 
of  the  u  variable,  but  no  resolution  of  the  variable  i. 


Example  4 .  Let  f(t)  =  J  s(t-nT), 

n  =  0 


where  s(t)  = 


1  for  t  e  [0,|] 


0  elsewhere 


f  is  a  coherent  train  of  N  + 1  pulses  of  width  — ,  and 

T 

separated  by  a  gap  -.  (Figure  3). 


Figure  3 


There  are  several  references  for  calculating  X .  The  most 
general  is  in  Bird  [2].  The  following  is  a  special  case 
of  Vakman's  calculation  [19], 

First  calculate  the  F.T.  of  f. 

F(v)  =  I  y  s(t-nT)e”^^^  dt  =  7  S(v)e^^^'^. 

J  n=0  n-0 

Then  from  (2.8)  we  see  that 


r  (u  ,v)  = 


F(v  +  ^)F(v  - 


- F-  ^  -i^T(m,+  n) 

I  I  S(v  +  .)S(v  -2)6  e 

n=0  m=0 


By  definition  (2.11),  the  AF  of  f  is 


X(u,t)  =  I  I 


^  ^  2iT  2  2 

n=0  m=0 


The  inner  integral  is  the  AF  of  s,  where  s  =  S, 


evaluated  at  t 


T(n-ir.)  i.e 


>;(u,T) 

Change  the 
k  =  n  -  rr; : 


X  (  u  ,T  ) 


H  M  -iyl(ni+n) 

I  I  e  X  (u,T-T(n-m) ) . 

n  =  0  in=  0  ^ 


variable  n  for  the  variable  k  by  letting 

then  k  takes  on  integer  values  to  M, 

’T  •  ukT 

r  /  V  2  r  -iir.uT 

=  I  X  (u,T-kT)e  Z  e 

k=-::  ^  m=o 


Define 


bk(u) 


'-‘I’"  fT 

^  V  -imuT 

I  e 

m=  0 


The  sum  is  easily  computed  so  that. 


-i~(k+N)  sin[u(N-t-l)T/2] 

fi,  (u)  =  e  .  uT 

sin  — 


and 


i: 

X(u,t)  =  I  X  (u,T-kT)<j!  (u)  . 

k=-N  ^  ^ 


Com.pare  the  mangitudes  of  ]({)  (u)|  and  |X  (u,0)|.  For 

K  S 

large  I.’,  |({),  (u)|  is  a  periodic  function  having  spikes 

of  magnitude  II  +  1 ,  (Figure  4), 


3  6.  Resolution  and  Arnbig’uity 

The  concepts  of  resolution  and  ambiguity  were  intro¬ 
duced  in  Section  3.  These  two  ideas  will  be  studied  in 
detail  in  this  section.  The  level  curves  of  the  AF  near 
the  origin  will  be  shown  to  always  be  ellipses  and  several 
examples  will  be  calculated. 


§6.1.  Resolution 

A  more  precise  definition  of  resolution  than  that 
introduced  in  Section  3  v;ill  be  used  here.  The  resolution 
of  f  is  the  vjidth  of  the  AF  at  the  origin  along  the 
coordinate  axes.  The  smaller  the  width,  the  better  the 
resolution . 

Since  we  are  interested  in  the  shape  of  X(u,t)  near 
the  origin,  let  us  approximate  X(u,t)  by  a  truncated 
2-dimensional  Taylor  expansion  about  (0,0).  Let  x  be 
an  arbitrary  AF  and  f  be  the  corresponding  generating 
function.  Let  subscripts  u,  and  t  denote  partial 
differentiation.  Then  by  the  definition  of  the  Taylor 
series  , 

(1)  x(u,T)  =  X(0,0)  +  X^_,(0,0)u  +  X^(0,0)t  +1  X^^(0,0)u^ 

+  X  (0,0)uT+^  X  (0,0)t^  +  ...  . 

^UT  2  TT 

It  is  assumed,  and  v;ill  be  verified  later,  that  for  the 

given  function  f,  X  =  X^  •  Equation  (6.1)  may  be 

normalized  by  dividing  through  by  x(0,0).  This  was  shown 

2 

in  (4.3)  to  be  I'f"  .  Henceforth,  assume  that  we  have 


r.crr.al  i  "eJ  the  function  f  by  '  f  ani  therefore 

\  ( : ,  i )  =  1 . 

Calculation  of  the  first  partial  derivatives  is 
strai rht forward  and, 

(2)  X,,(0,G)  =  -i  I  tl  f(t)|  ^  dt. 

Define 

(3)  a  =  Jt|f(t)l^dt 

then 

(4)  a  =  iX(0,0). 

u 

Ultinately  we  only  will  be  interested  in  the  magnitude 
of  X*  From  (4.6)  vie  know 


(5)  I  X  .(u,T  )  I  =  I  X.p  ^r(u,T  )  I  . 

a" 

Consequently,  without  loss  of  generality,  assume  the 
function  we  are  considering  in  (6.1)  is  f(t-a).  In  that 
case,  from  (4.6)  we  calculate  the  partial  derivative  of  the 
corresponding  AF  as. 


Similarly,  using  (2.11), 
(6)  X  (0,0)  = 

T 


dv . 


Define 


i  j  v|r(v)|^ 
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if  9 

(")  b  =  27f  j  ”  -iX^(0,0). 

Assume  that  the  F.T.  of  the  function  f  which  we  are 
considerinp  in  (6.1)  is  suitably  translated  by  b,  then 
it  can  be  shown  that, 

X  (0,0)  =  n. 

‘T 

Henceforth  assume  f  and  F  have  been  so  translated  and 
the  first  partial  derivatives  vanish.  (6.1)  may  now  be 
simplified  to 

(8)  X(u,T  )  =  1  +  I  x^^(0,0)u^  (0,0)uT  +  ^  (0,0)t^  +  .  .  .  . 

The  second  partial  derivatives  are 

(9)  =  ■{  t^lf(t)| ^  dt 

and 

(10)  x^^(0,0)  =  -  27  j  V^|F(V)| ^  dv. 

The  mixed  partial  derivatives  are 

(11)  x^^(o,o)  =  ~  j  v[F(v)rn:T7  -  F'(v)TrvT]dv 

and 

(12)  X^^(0,0)  =  Ij  t[f ’(t)f(tT  -  f(t)f ’ (t) ]dt . 

Notice  using  the  transform  pairs,  -itf(t)  •<— >■  F'(v)  and 
f'(t)  -»--»•  ivF(v)  and  Plancherel’s  theorem  that 

(13)  -ijv[F(v)FTvy-r’(vTrr7y]dv  =^jt[f '(t)f7TT-f(t)rTFT]dt. 


nence  X 


t  U  , '  ’ ) 


UT 


=  (0,0). 

T  'J 


Let  the  constants  in  (0.9),  (6.10)  and  (6.11)  be 
defined  -c^,  anl  -n  resrectivelv .  Then  the 

Taylor  expansion  up  to  the  quadratic  terr.s  in  (f.8)  r.av 
be  written. 


(14) 

X  (  U  ,T  ) 

= 

1 

19?  on 

|[d  u  +  2yut  +  D^-t  ^  ] . 

For  X 

equal  to 

a  constant,  this  is  the  equation 

ellipse 

in  the 

(  U  ,T  ) 

-plane.  That  is,  given  c  6 

the  level  curves 

of 

X 

near  the  origin  are 

(15) 

d^u 

^  +  2uut 

^2  2  2 
+  D  t  =  c  . 

It 

is  known 

that 

when  u  =  0  the  najor  and 

of  the 

ellipse  in  (6. 

15) 

lie  along  the  coordinate 

the  ( u 

, T ) -plane 

.  To 

see  what  conditions  need  to 

to  have 

U  =  0  , 

let 

f  ( t ) 

= 

and 

F(v) 

= 

A<v)ei'^'''h 

By  (6.12) 

y 

r 

=  -1 

ta 

^(t)6'(t)dt. 

by  (6.11), 


u  =  -  n—  vA^(v).f  ’  (  v)dv  . 


Hence  o 


0  if  f  is  real  or  F  is  real  . 
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We  showed,  in  Example  3,  a  signal  that  exists  for  all 
time  results  in  perfect  resolution  of  the  variable  in 
(doppler  shift).  This  is  equivalent  to  letting  the  pulse 
v;idth  of  Example  1  go  to  infinity.  The  consequences  are 
intuitively  satisfying  for  if  the  transmitted  signal  is  a 
single  frequency  for  all  timie , then  any  shift  in  frenuencv  o:" 
the  returned  signal  would  be  readily  apparent.  Conversely, 
an  attem.pt  to  resolve  the  time  delay  variable  from,  such  a 
signal  would  be  futile. 

The  last  exam.ple  concerning  resolution  will  dem.onstrate 
how’  to  sim.ultaneously  achieve  good  resolution  along  both 
axes.  The  idea  of  using  a  linear  FM  signal  to  achieve 
such  resolution  characteristics  was  a  major  development  in 
radar  technology.  Linear  FM  signals  solved  the  problem 
of  extending  the  range  of  a  radar  system  without  concurrently 
degrading  range  resolution. 


Exam.ple  5.  Let 


g(t) 


i-t^ 

9  ^ 

e  f  ( t )  , 


wnere 


f(t)  =  Prr,(t)  as  in  Examnle  1. 


Use  definition  (2.4)  to  calculate  the  kernel  function 


Y  ( t ,  T  )  as 
E 


laxt 


(16) 


Y  (t  ,t) 

*  IT  ' 


Y.(t  ,T  )  . 


hence  the  first  zero  along  the  i-axis  is  approximately 


3 


Consequently,  for  a  >  1  the  resolution  of  t  has  been 
irr.proved.  The  level  curve  for  given  c  and  large  T 
is  depicted  in  Figure  7. 


Figure  7 

§6.2.  Ambiguities 

The  AF  is  also  useful  in  characterizing  the  inherent 
ambiguities  in  a  given  radar  signal.  Should  a  given 
function  have  peaks  indistinguishable  from  the  peak  at 
the  origin,  then  identifying  the  point  which  corresponds 
to  the  actual  target  parameters  becomes  unclear. 

Consider  Example  4.  Here  f(t)  is  a  pulse  train. 

The  pulses  are  of  width  T/2  and  are  periodic  with 
period  T.  Intuitively,  one  might  reason  that  if  the 
first  pulse  leaves  the  receiver  at  t  =  0  and  does  not 
return  until  t  =  2T,  then  it  will  be  unclear  whether  the 


returned  pulse  is  a  reflection  of  the  first  pulse  trans¬ 
mitted  or  the  second  pulse  transmitted  at  t  =  T.  The 
reflected  energy  must  return  to  the  receiver  before  time 


t  =  T  for  the  receiver  to  unar.bi puous ly  r.easurc  the  ti 
delay  of  the  returned  pulse.  The  AF  (Firurc  5)  chara 
izes  this  uncertainty  because  it  has  peaks  alonr  the  t 
centered  at  t  =  FT,  k  =  . 

AF  ^Iso  characteri ses  uncertainty  in  the 
a  -  n .  'his  !act  is  net  so  reaiil'.’  arrarer.t  v.'hen  "7 
c;nri:c;  f(t)  in  isolation.  Hence  the  AF  rives  th 
v;:o.'_f_;-  ;eai^n^.r■  a  visual  tool  to  charact  er  i  se  ^  s  in.al  t  a 
e:  -£lv.  too  anoi  ruities  in  both  pararr.eters  . 


OD 

f(t)  =  I  5(t-nT). 

n  =  -a 


Figure  8 

To  derive  the  AF  for  f,  calculate  the  kernel  as 

OO  00 

Y(t,T)  =  [  6(t+|-nT)  I  6(t+^  -mT). 

n  =  -«>  m=-'» 

Y(t,T)  is  zero  except  when  t  =  kT ,  k  =  0,±1,...  . 


Therefore 


§  7 .  The  AT  as  a  Signal  Analysis  Tool 

The  AF  will  be  studied  as  a  general  analytic  tool 
in  this  section.  The  quadratic  phase  character  of  the  AF 
will  he  generalized  and  studied.  It  will  be  shown  hov;  the 
AF  characterizes  the  Heisenberg  uncertainty  principle  and 
several  examples  will  be  calculated. 

§7.1.  Decornposit ion  of  a  Linearly  Transformed  AF  [12]. 

Let  X  be  a  known  AF  of  a  given  function  f.  Let 
L  be  a  real-valued  matrix. 


Furthermore  let  L  map  the  coordinates  (u,t)  to  (U,T) , 

i .  e . 


Then  define 

(3)  X,(u,t)  =  X(U,T), 

i-j 

so  Xt^’-*’’'^)  "  X(au+bT  ,cu+dT )  . 

J-j 

Theorem.  Let  X  be  the  AF  of  f  and  define  Xj^  as  in 
(7.2).  X.  is  an  AF  if  and  only  if  det  L  =  1. 

Li 

proof .  First  the  necessity  of  the  condition  will  be 
verified  . 


Since  L(0)  =  0 


then  lx  (0,0)  I  =  |x(0,0)|.  Therefor 


from  (4.19)  we  conclude 


(4)  [ [ |x(u,t) I 

J  J 


dudT 


J  J 


X, ( u  ,T  )  I  dudT 


3y  definition  (7.3)  we  may  write 


(5)  II |X(u,T) I  dudT 
J  J 


J  J 


|X(U,T)  I  dUdT. 


Assume  the  Jacobian  does  not  vanish  in  IR^  then  we  may 
write 


(6) 


I X(u,T ) 1  dudx 


,2 

X(u,T ) I 


9(U,T) 
3  ( u  ,~t) 


dudT  , 


where  the  Jacobian  is 


3(U,T) 

aru  jTT 


det 


a  b 


c  d 


ad  -  be 


Thus  (7.6)  is  only  true  when  ad  -  be  =  1 . 

Next,  given  that  the  det  L  =  1,  it  will  be  shown 
that  X,  1  defined  in  (7.3)  must  be  an  AF.  This  will  be 

Lj 

demonstrated  by  decomposing  L  into  "elementary"  trans¬ 
forms  which  are  equivalent  to  operations  on  f  which  do 
not  alter  the  integrability  of  f.  Therefore  the  resulting 
function,  fj^,  may  be  used  to  generate  an  AF. 

Let  L  be  a  matrix  of  form  (7.1)  such  that  det  L  = 
Let  P  be  defined  for  some  a  €  IR , 


(7) 


P(a) 


Let  Q  be  defined  for  some  3  €  IR, 
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(S) 


Q(6) 


Then  if  c  0 ,  by  substitution  it  can  be  shown, 

(9)  L  =  P(a)Q(6)P(Y)  , 

where  a  =  ^(1-a),  B  =  -c  and  y  =  "'(l-d). 

If  b  0 ,  then  it  can  be  shown  that 


(10)  L  =  Q(a)P(6)Q(Y)  , 

where  a  =  ^(1-d),  B  =  -b  and  y  = 

If  both  b  =  0  and  c  =  0,  let  a  =  then 


(11) 


L(k) 


0 

k 


We  will  show  that  for  some  A  €  IR, 


(12) 


L(k) 


Q(a)P(B)0(Y)P(-A ) , 


where  a  =  r^( k-1),  6  -  kX  and  y  -  ,,(,  -1).  This 
kX 

be  shown  by  fixing  X  €  Pv  and  noting 


kX  k 


1 

k 


(13) 


L(k)P(X) 


-k;^ 

k 


This  is  in  the  form  for  the  decomposition  (7.10).  Eq 
(7.12)  follows  by  noting  that  P  ^(X)  =  r(-M. 


may 


uat icn 


Thus  L  can  be  decomposed  into  elementary  transforms 
of  the  form  P  and  Q.  Examiple  5  demonstrated  that 
transform  P  was  equivalent  to  m.ultiplying  f  by  a  linear 
FK  signal.  That  is,  given  a  function  f,  its  AF,  X, 
and  the  transformation  P(a)  as  defined  in  (7.7),  then 


(14) 
and 

(15) 


X  (u  ,T  )  =  X  ( u-ax  ,T  ) 

P 


q+  2 


f  (t) 

D 


f  (t)e 


ift 


Similarly,  given  f,  X  and  the  transform  Q(6) 
defined  in  (7.8),  we  will  show  that 


(16) 

and 

(17) 


Xq(u,t  ) 


X  (u,t-Bu)  , 


•  1+2 
i-rrt 


Use  the  F.T.  pair  e 


(2ti6)  ^  f(t)  *  e  . 


•  1^2  .3  2 

1 — "t  — i“V 

/i2TT"3  e  ^  and  the 


convolution  theorem  to  calculate  the  F.T.  of  f^  as 


Fq(v) 


(i)^  F(v)e  ^ 


Then  by  definition  (2.11),  the  AF  of  f  is 

q 


=  -4jT(v  +a)F(v-pe 


.3.  u>.2  .3,  .u,2 

-_lj(V  -  j  5) 


X ( u , T-u3 ) . 


are  still  L  -functions  therefore  X 


is  well  defined  and  the  theorem  is  proved. 

An  interesting  application  of  the  theor:  m  is  in 
property  7  .  There  we  saw  that  if 


L 


(18) 

then 

(19) 


X  J^(  U,T  ) 


f  (at) 


1  ,u  . 

T - [■  Xn  ^ 

|a|  fa 


This  is  transforjm  (7.11).  Therefore,  according  to  (7.12), 


.X  .  2  .1.2  .6^2  .1^2 

-i-t  ift  1— t 

(20)  fiat)  =  A(((f(t)e  ^  )*e  )e  ^  )*e  , 


where  a  =  (a-1),  B  =  aX  ,  y  =  ^  is 

aX  aX  a 

arbitrary,  and  A  is  the  constant  of  proportionality 
dependent  on  a  and  X . 

Example  7 .  From  (2.13)  we  know  Xp(u,T)  =  2it  X^(t,-u). 
Thus  the  AF  of  F  results  in  the  linear  transform  of 
the  AF  of  f  represented  by. 


L 


1 

0 


From.  (7.9)  we  may  decompose  L  into, 

L  =  P(-1)Q(1)P(-1) . 


Therefore  a  real-time  spectrum  analyzer  may  be  constructed 
as  follows; 


1 


w 


c 


Figure  10 


Central  Moments  and  the  Uncertainty  Principle 


Signal  duration  and  bandwidth  are  frequently  used 
quantities  in  signal  analysis.  Perhaps  the  most  common 
way  to  characterize  these  features  is  in  the  moments  of 
inertia.  In  the  following,  the  central  mom.ents  v;ill  be 
shovjn  to  be  the  derivatives  of  the  AF  at  the  origin. 

The  shape  of  the  aF  will  be  related  to  the  Heisenberg 
uncertainty  principle  and  several  interesting  examples 
will  be  presented. 

In  most  signal  analysis  applications,  f(t)  is  the 

complex  representation  of  the  real-valued  signal  of 

interest.  Therefore  the  central  moments  will  be  defined  in 
2 

term.s  of  1  f I  .  Define  the  energy  of  the  function  as 


(21) 


I  I f (t ) 1 ^  dt . 


3y  definition  (1.3)  and  equation  (4.3)  we  also  note 


X  (  0 ,0)  . 


The  center  of  gravity  is  defined 


■  t|f(t)I^  ^ 
'  ‘n 


V.'e  vjill  henceforth  assume  mi^  =  0.  This  is  a  valid 
assumption  because  we  may  always  translate  f  to  make 


=  0 .  Finally,  define  the  moment  of  inertia  or  effective 


signal  duration  as. 


This  was  defined  as  d  in  Section  5.  Pv  definition 
(1.3)  we  may  also  write 

(25)  =  — 'ltf(t)''^  =  d^. 

.  mg 

Similarly,  the  moments  of  F  =  f  are  defined; 

(26)  Mg  =  ~  j|F(v)l^  dv  =  II  F' ^  =  X(0,0), 

(27)  ~  ^  J  v|F(v)|^  =  0  (by  assumption), 

and 

(28)  M^  =  ~  I  v2iF(v)|2  =  1  !!vF(v)II^  = 

'  j  0  0 

Recalling  equations  (6.9)  and  (6.10),  we  see  that 

(29)  m-  =  -X  (0,0) 

2  uu 

and 

(30)  M^  =  -X  (0,0) 

2  TT 

Thus  the  width  of  the  AF  function  along  the  x-axis 
depends  on  the  spectrum  of  the  generating  function  and  the 
width  along  the  u-axis  depends  on  the  duration  of  the 
signal.  The  Heisenberg  uncertainty  principle  (UP)  is 

(31)  i'tf(t)  :vF(v)  li^  g 


Th-3  the  curvature  of  x  at  t];e  oririn  is  constrained 
bv  the  i:F; 

( ?  ^ )  X  C  ,  ■  ) X  ( C  ,  G)  >  i 

UU  TT  / 

Unlike  the  Wifner  distribution  of  the  next  chapter,  this 
does  not  restrict  X  frorr  becoming  concentrated  about 
the  oririn.  Klauder  demonstrated  a  compressible  AF  in 
1960  [7]. 

We  next  consider  the  signal  duration  and  bandwidth 
of  the  function  f  .  f.  is  the  generating  function  of 

L  lj 

the  linearly  transform.ed  AF,  Xj^,  defined  in  (7.3).  To 
simplify  notation,  given  f,  let  the  moments  of  inertia 
of  f  and  F  be 


(34)  vj  =  m^  =  (from  (7.29)), 

and 


(35)  W  =  =  -X  (0,0)  ( from  (7.30)), 

I  TT 


respectively . 

Recall  in  Section  6  we  defined 


y  =  X  (0,0)  =  X  (0,0)  . 

UT  TU 


and 


tiate  (7.3)  to  shovi 


(  36) 


w. 


2 

a  w 


+  2ac  u  + 


c^W, 


(37)  =  b^w  +  2bd  u  +  d^W, 
and 

(38)  y  =  abw  +  (ab+bc)y  +  cdW. 

Li 


These  relationships  will  be  applied  in  the  following 
exanples . 


2 

Example _ 8.  Let  r  €  L  (IR)  and  real-valued.  Define 

i-t^ 

(39)  f(t)  =  r(t)e  ^  . 

The  F.T.  of  f  is 


(40)  F(v) 


r(t  )e 


r(t  )e 


This  is  not,  in  general,  computable  in  closed  form. 
Usually  a  is  assumed  large  enough  so  that  stationary 
phase  arguments  may  be  applied.  Then 


(41) 


F(v) 


1  , 
^2a 


r(-)  . 


This  may  be  used  in  (7.28)  to  approximate  the  effective 
bandwidth  of  F  as 


(42)  VJ 


1  f 


^  v"|F(v)rdv  =  - 

F  2it  '  '  a 

J  < 


2  I  /  ^  \  I  ^  j  ^ 

v  r  ( - )  dv  =  a  w  , 
'  n  '  r 


where  w 


was  defined  in  (7.24). 


>^i< 


We  may  use  the  results  of  the  previous  analysis  to 
get  an  exact  expression  for  effective  bandwidth.  From 
equation  (7,15)  and  (7.7)  we  see  that  the  transform 


matrix  is 

P(a)  , 

and  a  =  1,  b  =  -a, 

Hence  from 

(7.37) 

we  see  that 

(43) 

2 

=  a  w  -  2au  +  W 
r  r  r 

r  is  real- 

valued 

,  hence  u  =0  and 
r 

(44) 

“l 

a  vj  +  W  . 
r  r 

For  large  a,  the  equation  (7.44)  agrees  with  the  approxi¬ 
mation  (7.42) . 

Consider  equation  (7.43)  as  a  quadratic  in  a.  For 

fixed  r  (not  necessarily  real-valued)  w  ,  y  and  W 

r  r  r 

are  fixed  and  (7.43)  is  a  parabola  opening  upward. 

For 


u  =  0 ; 

min  W^(a)  = 

VJ  at  a  = 

0  ,  and  for 

r 

F 

r 

y  0 ; 

min  W„(a)  = 

W  -  y^/w 

at  a  =  y  /w  . 

r 

F 

IT  IT* 

r  r 

Hence  we  see  the  effective  bandwidth  decreases  for  a 

U  >0 
r 


( Figure  11 )  . 


Figure  11 


Exar.ple  9.  Let  r  ^  L^(IR)  and  real-valued.  Define 


(45) 


f^(t) 


•  1+2 
loot 


[r(t)e*2^  3  - 


Using  (7.7)  and  (7.8)  v/e  see  that 


(46) 


P(-a)Q(3) 


l-aB  a 
-S  1 


and 


(47)  X.(u,t)  =  X  (u-aSu+ax ,-Bu+t )  . 

Lj  ip 

r  is  real,  so  u  =  0.  Fix  a  ^  IR  then  by  (7.36), 

(48)  w  =  Wt  ( B )  =  (1-aB)^  w  +  B^  17  • 

L  L  r  r 


Expjanding  we  obtain 


(49) 


(a  +V.’  )2  -  -  2a V7  0  +  vj  . 
r"  r-  r  r 


This  is  a  quadratic  opening  upward  with  a  ininimur.  at 


(  50) 


aw 

B  .  = 

'"i"  tv; 

r  r 


For  a  >>  ,  this  may  be  approximated 


as 


(51) 


min 


_ 1/a _ 

1+W^/a^w^ 


1 

a 


Hence,  we  see  that  although  a  signal  disperses  when  passed 
through  a  quadratic  phase  filter,  in  this  case  we  get  some 


signal  compression  (Figure  12). 


Figure  12 


From  (7.51)  and  (7.i-8)  we  find  the  minimum  sitnal  duration 


III  V’lGI.'ER  DISTRIBUTIOri 


In  this  chapter  the  VJipner  distribution  v.’ili  le 
presented.  It  will  be  defined,  its  salient  properties  will 
be  stated  and  proved,  and  exanples ,  which  '..’ill  serve  to 
illurr.inate  those  properties,  v;ill  be  calsuJated.  rinall-; 
sone  interestinp  results  concerninr  the  rorer.ts  or'  the 
Wiener  distribution  (WD)  will  be  presented. 

The  Wigner  distribution  was  originally  introduced  in 
the  context  of  quantum  mechanics  in  19  32.  It  v;as  reintro¬ 
duced  by  Ville  in  1948  as  a  tool  for  signal  analysis.  The 
r.ost  recent  appl •' cations  of  the  Wigner  distribution  have 
been  in  the  field  of  optics.  In  1983  Claasen  and  Mechlen- 
brauker  again  studied  the  WD  v;ith  regard  to  its  potential 
as  a  signal  analysis  tool  [3]. 

§  8 .  Def inition 

The  WD  is  defined  for  continuous  functions  f, 

mapping  the  real  line  to  the  complex  plane,  which  are  of 

2 

inite  energy.  If  f  ^  L  (R)  or  if  f  is  not  continuous 
then  f  -will  be  considered  a  peneralized  function  and 
operations  on  f  will  be  in  the  distribution  sense. 

The  WD  is  defined  as 

!  T  TT  -ivT 

(1)  W(t,v)  =  f (t  +  y)f (t  -  j)e  dt. 

I.'ote  that  the  definition  is  stated  as  an  "auto-Wigner 
distribution."  This  definrtion  easily  generalizes  to  a 
"cross-Wigner  distribution . " 


As  in  the  case  of  the  AF  we  define 


(2)  Y  ^t,T  )  =  f(t  +^)f(t  -^)  . 

Then  for  fixed  t,  the  WD  is  the  Fourier  transform  with 
respect  to  the  variable  t  of  the  function  Y(t,T),  That 
is  V  t  €  R , 

(3)  W(t,v)  =  Y  ("t  dT  . 

Let  F  =  f.  Then  it  can  be  shown  that 

(4)  W^(t,v)  =  ~  W^(v,-t). 

If  we  define 

(5)  r(u,v)  =  F(v  +y)F(v  -^) 

and  recall  (2.6),  which  defined  r(u,v)  as  the  two- 
dimensional  Fourier  transform  of  Y(t,T),  then 

(6)  Y(t,T)  h^W(t,v)  r(u,v). 

This  says  that  the  WD  is  the  intermediate  result  of  the 
iterated  two-dim.ensional  F.T.  of  Y(t,T)  when  the 
transformation  with  respect  to  t  is  taken  first.  It 
was  shown  that  for  fixed  v  €  IR  r(*,v)  €  L^(IR)  hence, 

(7)  W(t,v)  =  I  F(v  +-y)F(v  -  du 
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§9.  Properties  of  WB 

The  properties  will  be  listed  in  somewhat  the  same 
order  as  in  the  previous  chapter.  Iv’here  proof  of  the 
property  is  similar  to  the  arpument  for  the  AF ,  the 
property  vjill  be  stated  without  proof. 

1.  Global  Maximum 

The  WD  has  global  extrema  at  the  origin  only  if  f 
is  real  and  sym.metric.  This  property  does  not  necessarily 
generalize  to  the  "cross-Wigner  distribution."  If 
f  :  1?. -»■  P.  and  f(t)  =  ±f(-t)  then 

(1)  |W(t,v) [  S  |W(0,0)  1  . 

This  property  is  stated  in  terms  of  the  modulus  of 
the  WD.  Although  the  WD  is  always  real  (property  2), 
its  value  at  the  origin  may  be  negative,  e.g.,  the  case 
when  f  is  odd. 

2 .  Symmetry 

a.)  The  V/D,  regardless  of  the  nature  of  the  function 


f ,  is 

always  real  valued. 

That  is 

(2) 

W( t ,v)  = 

W( t , v) . 

Proof . 

Note  that  y’(  t ,  t7  = 

Y(t,-T).  Then  by  definition 

(8.3), 

=  YTt,r)e  dx  = 

. 

1  y(t  •>“T  dt  =  W(t,v). 

This  property  generalizes  to  the  "cross-Wigner  distribution 


since 

b.)  If  f  is  real-valued  and  symmetric  about  the 
origin  then, 

(3)  W(-t,-v)  =  W(t,v). 

Proof .  Consider  only  the  case  f(t)  =  -f(-t).  Then 

Y(_t,T)  =  f(-t  +|)f(-t  -^)  =  f{t-|)f(t+|)  =  Y(t,T). 

The  property  follows  from  definition  (7.3). 

3.  Translation 

a.)  Translation  of  the  function  f  results  in  a 
translation  of  the  WD.  Let  T  be  the  translation 
operator.  Then 


(4) 


f(t,v)  =  VJ^(t-a,v) 
a 


Proof .  Notice  that  ^(t,T)  =  y^Ct-a,!).  Consequently, 


^  ( t ,  V )  = 
a 


y^(t-a  ,v)e  dr  =  W^(t  -a,v). 


b.)  Translation  of  the 
a  translation  of  the  WD  of 

T  F(v)  = 
a 

(5)  and  h(t)  *-*- 


F.T.  of 
f .  If 

F( v-a) 

F( v-a)  , 


f  also  results  in 


then 


W  ( t ,  V ) 


W^(t ,v-a)  . 
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Proof.  It  is  known  that  h(t)  =  e  f(t).  Hence 


Wj^(t,v)  =  j  f(t  +^)f(t 


^)f(t  -^)e 


ia(  t  +  i  )  -ia(  t  -  1) 


2 

e  dT 


=  j  f(t  +^)f(t  dx  =  W^(t,v-a) 


4.  Modulation  by  a  simple  sinusoid. 

The  previous  property  may  also  be  stated  in  terms  of  a 
modulation  operator.  That  is,  let 


M  f(t)  =  f(t) 

a 


Then  (9.5)  may  be  written 


3. 


VJ^(t ,  v-a) 


Similarly  if  M  F(v)  =  e~  F(v)  and  h(t)  -«->■  e' 

"■  5l 

then  (9.4)  may  be  written 


W^(t,v)  =  W^(t-a,v). 


F(v) 


An  interesting  computational  rule  is  a  result  of  (9.4) 
and  (9.6).  Combining  these  two  properties  we  get 


b  a 


W^( t-a , v-b) . 


Let  t  =  V  =  0  and  change  the  sign  of  the  parameters  a 
and  b .  We  obtain 


T  f  ^  0  >  0 ) 
-b  -a 


W^(a,b) . 


r  . 


Hence,  the  wD  of  a  function  f  may  be  evaluated  at  any 
point  (a,b)  in  the  plane  by  first  translating  and  modulat¬ 
ing  f,  then  evaluating  the  WD  of  the  resulting  function 
at  the  point  (0,0). 

5.  Multiplication  of  two  functions. 

If  h(t)  =  g(t)f(t)  then 

•)  f 

(10)  W  (t,v)  =  (t ,r )W^( t ,v-r )dr . 

h  2it  J  g  ’  f  ’ 

0.  Filtered  Functions. 

Let  h(t)  =  (g*f)(t),  where  g  is  the  impulse  response 
to  a  linear,  time-invariant  filter.  Then 

(11)  W,  (t,v)  =  [  W  (r ,v)W^(t-r ,v)dr . 

h  ’  1  g  ’  f 


7.  Invertibility  and  Uniqueness. 

From  (8.3)  it  is  apparent  that  given  a  WD  W,  one  can 
uniquely  invert  the  F.T.  to  recover  the  kernel  function 
y.  Hence,  given  W(t,v)  is  a  WD,  then 

(12)  f(t  +^)f(t  -j)  =  —■  J  W(t,v)e^'^^  dv. 

Let  t^  I  ’  ^2  ~  ^  ~  \  that 


1  f  t.+to  iv(t  -t  ) 

_  I  dv. 


If  t  =  t^,  t^  =  0  then 


,,,t  ,  ivt 

y(-,v)e 


(14) 


f(t)TTTy 


lienee  the  generating  function  f(t)  can  be  recovered  to 


within  a  constant  f(0).  As  in  the  case  of  the  AT,  it 
can  be  shown  that  this  constant  ir.ust  have  modulus  1 . 
Furthermore,  by  letting  t^  =  t^  =  t  in  (9.13)  we  have 
the  unusual  result  that  Vt  k  IF 


(15) 


lf(t)l^ 


-  -  I  V.’  ( t  ,  V )  d  V  . 

2ti  J 


The  novelty  of  (9.15)  lies  in  the  fact  that  the  WD  is  not 
always  positive.  Yet  (9.15)  says  that  regardless  of  the 
nature  of  f,  at  any  particular  t,  the  integral  of  the 
WD  over  all  frequencies  results  in  a  nonnegative  number. 


8.  Volume  Invariance. 

(16)  “  W(t,v)dtdv  =  I'f' 

/TT 

2 

where  the  L  -norm  of  f  was  defined  in  equation  (1.3). 
This  is  a  direct  result  of  (9.15). 


9.  Analytic  Signals, 

Let  f  be  the  complex  representation  of  the  real 
signal  s(t)  as  in  (A. 3).  Then 


(17)  W^(t,v) 


V  -  0 

V  <  r , 


where  h(r,v)  =  ^(sin  2vr)  .  Equation  (^.17)  is  th.e 
convolution  of  the  WD  of  the  real-valued  signal  s(t) 
with  the  Fourier  kernel  h(r,v). 


Proof . 


Recall 


( A . 6 )  to  see  that 


2S(v+^)  for  u>  -2v 


(18) 


r(v 


S(0) 

0 


u  =  -  2v 
u  <  -  2  V 


and 


2S(v  -  ~)  for  u<  2v 


F(v  -^) 


SCO) 

0 


u  =  2v 
u  >  2v . 


Temporarily  assume  v  >  0.  Hence,  using  definition  (8, 
we  see 


(19)  W^(t,v) 


2  .u 


-2v 


S(v+^)S(v  -  |)e^^^  du, 


Let 


P2v‘“> 


1  for  |u  <  2v 


0  elsewhere. 


Then  (9.19) is 


(20)  h'  (t,v) 


^  f  -  (u)S(v.-pS(v-H)el“P 


2v 


du , 


It  is  known  (— -)  sin  2vt  p„  (u).  By  the  convolution 
TTt  t  ^2v 


theorem,  since  W  -*->■  f ,  we  conclude 

s  u 


(21)  W^(t,v) 


4 

tt' 


(F )  sin(  2vr  )VJ  (t-r,v)dr,  v  >  0. 
r  s 


Should  V  <  0  then  (9.18)  states  that 


F(v  +^-) 


0 


for  u  <  -2v 


F(v  -^)  =  0  for  u  >  2v. 

Hence  F( v  +  ^) F( v  -  =  0  Vu  ^  IR,  and  the  property  is 

verified.  * 

§  10 .  Examples 

Example  10.  Let  f(t)  =  p^(t).  Then  refer  to  Figure  1 


and  note 

that 

the  limits  of  integration  of 

the 

WD 

for 

fixed  t 

are 

T  =  -2T  +  2|tl  to  T  =  2T 

-  2|  t| 

• 

Then 

for  T  > 

0, 

f2(T-l tl )  .  „  f (^)sin 

2v(T-l 

t]  ) 

for 

1  tl  <T 

(1)  W(t 

,v)  = 

J-2(T-|t|)  Q 

for 

|tl>T 

Note  that  since  f  is  real  and  even,  the  WD  has  a  global 
maximum  at  the  origin.  Except  for  a  scale  factor,  the  -WD 


is  the  same  as  the  AF. 
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Exar.Dle  11. 


Let 


f(t) 


I  t|  <  T 


elsewhere . 


From  (9.10)  and  (10.1)  we  conclude 


(2)  W(t,v)  =  { 


( v-a ) 


sin[  2  ( v-a )  (T-l  tl  )  ]  for  1  tl  <  T 


elsewhere . 


Example  12 . 
Let  f(t) 

Note  that 


Ae^^^,  V  t  €  IR  and  some  A  €  C. 
Y  (t  ,T  )  =  I  A|  ^  . 


Then  formally, 

(3)  W(t,v)  =  I  A|  ^  J  dx  =  2tt  |  Aj  ^  6  (v-a), 

v;here  6(v)  has  been  previously  defined  in  (1.7). 


Exar.p  • 

Let 


.  a  2 
i  «t 

f^(t)  =  e  , 


.  a .  2 


t  <  T. 


Define  g(t)  =  e  ^  'and  f(t)  =  p^(t).  We  will  find  the 
Vj’D  for  r(t)  for  all  time,  then  use  Example  10  and 
property  5  to  find  the  WD  of  f ^ .  Formally, 

•  3.  f  .  t,2 


v;  (t,v) 
r 


r  i‘^(t+i)"  -i?(t-i) 


2 '  "  ■  2 '  '"2  '  2  ■  -IVT  , 

e  e  e  dx 


J 


latx  -ivx 


dx  =  211 6  ( v-at )  . 


Recalling  (10.1)  and  (9.10)  we  conclude 


l\  (t,v) 
^1 


---  J  271(5  (r-at)^----^  sin  [  2  ( v-r )  ( t-|  T|  )  ]dr 

{  ( — )  sin[  2(  v-at )  (t-|  Tl  ]  I  t|  <  T 
I  v-at  '  '  II- 


elsevihere  . 


This  may  alsc  be  written. 


VJ  (t,v) 
^1 


W^( t , v-at ) . 


We  see  that  multiplying  a  function  by  a  linear  FM  signal 
results  in  a  linear  transform  of  the  original  V.’D. 


Example  14 .  Let  f  be  given,  and  W(t,v)  be  the  VJD  of 


f.  Define 


f  (t) 

1-j 


--■  i-^t^ 

(2TTb)  2  f(t)  fte 


V.e  will  calculate  the  WD  of  f  . 

Lj 

.  1  ^2 

Using  the  transform  pair  e 


/2"tt  ib  e 


-i^v 


we  can  formally  calculate  the  FT  of  f  as 


Fl(v) 


1  .b  2 

2  ■^2'' 

(i)  e  F(v). 


Use  the  definition  (8.7)  and 


wpt,v)  --  A  j  rpv*|)F~(v.|Tel“"  du 

=  ril  rCvdbnTll)  du 

2itJ  2  2 


( t-bv ,t )  . 


§  11 . 


C entr al  Honents 

It  is  knov;n  that  the  central  rricr.ents  of  a  function 
provide  sone  insipht  into  the  shape  and  character  of  a 
function  without  ^ully  des  ur :  r  inp  it.  The  follov;ing  v;ill 
generalize  the  notion  o:  ccr.*r:il  Tr:oTrents  in  IP.  to  local 
and  global  mon.ents  in  IP  ‘  .  V.v  'will  see  that  these 
ir.oir.ents  of  the  WD  are  sor.e  v.’ell  known  quantities  in 
signal  and  netvjork  analysis. 


Let 

the  central  mom.ents  of 

|f|^ 

be  m.,.^  ,  m^ 

and 

r.j  ■ss 

defined  in 

(7.21),  (7.2?) 

and 

(7.24).  Let 

■  the 

moments 

of  |F|^ 

be  ,  K,  and 

.1  1 

w 

"2 

as  defined 

in  (7.25) 

(7.28). 

We  will  first  study  the 

local 

momients  of 

the  WD . 

For 

’  fixed  t. 

the  local  ave 

rape 

of  the  WD 

with 

respect 

to  V  has 

been  calculated  in 

(9.15)  as 

(1) 

nyt)  = 

r 

j  VJ(t,v)dv  = 

1  f(t) 

if 

This  is 

the  power 

in  f  at  time 

t . 

The  local  center  of 

gravity 

with  respect  to  v  is 

(2) 

npt,  =  1 

vW(t,v) 

This  may  be  formally  calculated  using  Fubini's  theorem  and 
the  F.T.  pairs  v  27t6'(t)  as 
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(3) 

n^(t)  =  J  [f(t+|-)f(t-^)]6’(t) 

=  ||[f '(t +^2^f(t -l)-f(t +I)f' (t -^)]5:t)  n-j7-|y^ 

=  277^ t7  tf 'Ct)fTt7-f(t)f’rt7]. 

If  f  is  real-valued  then  n^(t)  to.  If  f  is  complex- 
valued  then  (10.3)  can  be  put  in  a  more  meaningful  form, 
notice  that  since 


[ff ’  -  ff ’  ]  =  Im  ff ' , 

then  (10.3)  is 

(U)  n^(t)  =  In'.  =  Im{-^-  In  f'(t)}. 

Therefore  if  f(t)  =  a(t)e^*^^'*'^  we  see  that 
(5)  n^(t)  =  0'(t). 


For  complex  representations  of  a  real-valued  signal 
(Appendix  A)  this  is  the  instantaneous  frequency  of  the 
function  f.  The  VJD  permits  the  generalization  of 
instantaneous  frequency  to  arbitrary  complex- valued  function 
The  local  moment  of  inertia  for  fixed  t  is 


(  6  )  n^ ( t ) 


±  f 

2tt 


r 


( v-n^ ( t ) )  W( t , v ) 


A  j  K(t,v) 


_dv__ 

n  ( t ) 

0 

n^(t)^ 
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The  WD  is  not  always  positive  so  interpreting  this  as 

the  spread  of  the  VJD  is  not  entirely  accurate.  To 

derive  some  meaning  from  n^Ct)  we  must  put  (10.6)  in  a 

2 

different  form.  Use  the  F.T.  v  2tt5"(t)  to  formally 

calculate  n^d)  in  terms  of  f.  We  see  that 

(7)  n^Ct)  =  -~^py[f"(t)fTtT-2]  f ’(t)|  ^  +  f(t)F'dl]  -  n^(t)^. 

0 

It  can  be  shown  that 


(8) 


n^Ct) 


1  d  f’(t), 
'  2  ^^^dt  f(tl-^ 


If  f(t)  =  then 


(9) 


n^Ct) 


dt ' 


\  a('t)  I 


Hence  the  local  second  moment  of  the  WD  is  independent  of 
the  phase  of  f.  Furthermore,  for  any  a  6  IR  and  any 
C  >  0, 

(10)  n2(t)  =  0  if  and  only  if  a(t)  =  . 

The  sufficiency  of  the  condition  is  shown  by  substitution. 
The  necessity  of  the  condition  is  shown  bv  letting 

,2 

(11)  inla(t) i  =  0 , 

dt 

and  integrating  twice. 

Consider  now  the  local  moments  of  the  V.’D  with  resmect 
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Fix  V,  then 


r 

(12  ) 

1 

< 

M 

I  W(t,v)dt  = 

1 F(v) 1 

3_.':;ildrly ,  the  local  center  of  gravity  is 
(i:)  =  |  tv.'(t,v) 

Tnis  can  be  shov;n  to  be 

(14)  ’i^Cv)  r  2^Ttl  (v) ] 

(15)  h'^(v)  =  =  -Im  {^lnF(  v)  }  , 

If  F(v)  =  A(v)e'^^'^^  then 

(16)  N^(v)  =  <f'(v). 

Should  F(v)  be  a  syster.s  transfer  function,  then  (11.16) 
is  the  group  delay.  Again  the  V,’D  perrdts  a  specific 
definition  to  be  generalized. 

By  similar  calculations,  as  in  deriving  (11.9)  it 
can  be  shovjn  that  the  local  moment  of  inertia  with  respect 
to  t  is 

(17)  hd(v)  =  -  --  — ^lr:|A(v)|. 

dv 

The  global  moments  of  the  VJD  are  taken  over  the 
entire  (t,v)-plane,  and  will  be  shown  to  be  moments  of 
the  generating  function,  f,  and  its  F.T.,  F.  The 


global  average  of  the  WD  was  previously  calculated  in 


(9.16)  as 


(18) 


“  J I  1  v)dtdv 


Similarly 


(19)  Nq  =  2^  J J  W(t,v)dtdv 

Hence  n^^  =  . 


n  f" 


2 


T'': 


2 


The  global  mean  or  center  of  gravity  with  respect  to  the 
variable  v  is 


(20)  n 

1 


_1_ 

2Tr 


’ 

. 


vW ( t , v ) 


dvdt 


v| F(v) 


2 


dv 

2  ■ 

I!  F!i 


The  last  equality  results  from  interchanging  the  order  of 
integration  and  (11.12).  Hence  the  global  center  of 
gravity  of  the  VJD  with  respect  to  the  frequency  variable 
V  is  the  same  as  the  center  of  gravity  of  |F(v)|  ,  where 

F  is  the  F.T.  of  f. 

In  a  similar  manner,  the  global  mean  of  the  WD  with 
respect  to  the  variable  t  is  defined  as 


(21)  i;. 


JJ 


‘  t|f(t)|^ 


N 


0 


The  global  moments  of  inertia  are  defined  as 


(22) 


f  f 


f  --  ^2,,..  .dvdt 
( v-n  )  W(t ,v)-- — 


|(  v-n^  ) 


F(v) 
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and 

(23)  ~  I  I  (t-lT  )^W(t,v)-^--  =  I  (t-fl,  I  f(t)  I  ^ 

2  27:  jj  1  j-  J  1  '  '  ,2 

'o  ■  ■ 

We  see  that  the  moments  of  inertia  are  non-negative  and 
hence  may  be  accurately  interpreted  as  a  measure  of  the 
spread  of  the  'WD. 

§  1 2  The  Uncert a i nt y _  Fr incipl e 

The  Heisenberg  uncertainty  principle  (UP)  constrains 
the  mom.ents  of  inertia  of  a  function,  f,  and  its  F.T.,  F. 
Consequently  the  U?  must  also  constrain  the  global 
moments  of  the  WD.  A  form  of  the  UP  is 

(1)  'ltf(t)il!!vF(v)ll  >  |!;fll^. 

In  the  case  of  non-centered  moments  (12.1)  may  be 
also  written 

(2)  ■(t-a)f(t)l!!'(v-b)F(v)!!  >  |l!fl|^, 

where  a  and  b  are  the  respective  centers  of  gravity. 
Using  similar  notation,  the  global  moments  (11.22)  and 
(11.23)  may  be  written 

(3)  =  IKv-H^  )F(v)!'^  — N  =  I' ( t-N  )  f  ( t )  !i  ^ 

/  X  j!  p  !j  ^  i.  J-  |i  ^ 

Without  loss  of  generality,  assume  f  and  F  have  hence¬ 
forth  been  shifted  so  that  n^  =  =  0  .  VJe  will  now 

show  that  the  WD  cannot  be  concentrated  arbitrarily  close 


to  the  origin.  Let  x  and  y  be  fixed.  Then  it  can  be 
shown  that , 

2  1 

(4)  ir.inCp  x  +  — 2  y)  =  2/3?/. 

P  P 

.2  ,2 

If  X  =  j,tf(t).’  and  y  =  !'vF(v)‘'  then  this  irinlies  that 
for  any  p, 

(5)  p^:!tf(t):;^  +  -^:vF(v)r^  >  2::tf(t)^:"  vF(v)  .. 

p 

Then  by  (12.1)  we  conclude  that 

(6)  p^lltf(t)|!^  +  -^'■'vF(v)!l  ^  >  llfl'^. 

P 

Using  definitions  (12.3)  with  n^  =  =  0  this  may  also 

be  written 

(7)  p^TT.  +  ii  >  1. 

P 


In  terms  of  the  WD,  this  is 


~  JJ(p^t^  + -^  v^)W(t  ,v)dtdv  >  JJ  V/(t,v)dtdv. 


(12.8)  may  have  the  following  interpretation:  For  fixed 
p  ^  P,  let 


(1)  g(t,v)  =  p^t^  +  . 

P 

2 

Then  Vr  e  P,  g(t,v)  =  r  is  a  weight  function  in  the 

2 

integral  (12.8),  v;hich  assigns  the  value  r  to  all 


values  of  the  WD  which  lie  on  the  ellipse. 


Hence  f[(t,v)  suppresses  the  values  of  the  VJD  near  the 
origin  and  amplifies  the  contribution  of  the  WD  away  from. 


the  origin.  The  inequality  in  (12.8)  means  that  the  V.’D 

cannot  be  totally  concentrated  in  an  arbitrarily  small 

region  about  the  origin.  If,  for  example,  the  VJD  were 

2 

to  vanish  off  an  ellipse  such  that  r  <1,  then  inequality 
(12.8)  would  not  hold. 

The  inequality  (12.8)  does  not  preclude  a  highly 
concentrated  VJD  which  has  a  small  contribution  far  from 
the  origin.  If  the  WD  is  to  characterize  the  Heisenberg 
uncertainty  principle,  then  we  must  also  prohibit  just 
such  a  WD.  To  see  that  such  WD’s  are  in  fact  impossible 
we  need  the  following  theorem  [5]. 

If  a>0,b>0  and 

TT.  .27Tx  v2 

,  f  f  — (t-r)  --(v-s) 

(11)  W^b(t,v)  =  — ^  ^  W(r,s)drds 

2^ / ab  J  J 

is  the  Weierstrass  transform  of  the  WD  W(t,v),  then 

9  r  2 

(12)  Vt,v  ^  IT,  W  ,  (t,v)  s - ^ —  |f(t)|  dt. 

1+2/ab  J 

The  proof  of  this  statement  depends  upon  expanding  f  in 

2 

an  orthogonal  system  in  L  CT)  which  is  related  to  Herm.ite 
Dolynomials  and  is  beyond  the  scope  of  this  paper.  VJe  mav 
use  this  result  for  the  specific  values  t  =  2r  and  v  =  2s, 
then  (12.11)  is 


e 


e 


W(r , s ) drds  . 


(13)  W  ,  (2r,2s)  = 

an 


2tt  /ab 


(12.12)  then  becomes 


r  r 


2  2 
TT  r  TT  S 


(14)  A  jj  e  ^ 


^  VJ(r,s)drGs  5 


f(t)|  dt 


1+2/ab  J 


Notice  that 


(15) 


__2/ab_ 

1+2/ab 


(  1 - ) 

\  1  +  2/ab  ' 


hence  (12.14)  is 

(16)  jj(l 

VJe  now  have  an  inequality  using  a  weight  function  which, 
for  fixed  a  and  b,  goes  to  1  with  increasing  values 
of  r  and  s .  This  precludes  a  WD  as  was  previously 
described.  Therefore  from  (12.8)  and  (12.16)  we  conclude 
that  the  WD  cannot  be  arbitrarily  concentrated  about 
the  origin  (or  in  the  case  n^  ^  0  or  ^  0 ,  about 

the  center  of  gravity). 


-  e 


.IE_ 

a 


ITS 


'W  ( r ,  s )  dr  d  s 


1+2/ab 


VJ(  t ,  v)dtdv 


IV.  EPILOGUE 


In  the  previous  chapters  two  time-frequency  functions 
were  studied.  Their  properties  were  reviewed  and  some 
applications  shown.  Attention  was  directed  to  their 
relationship  to  the  Heisenberg  uncertainty  principle. 

In  the  next  section  we  will  state  the  relationship  between 
these  two  transforms  and  consolidate  their  similarities  and 
differences . 

S 1 3  .  up _ ^vs_^ _ ^ 

Recall  equations ( 2 . 9 )  and  (8.6)  and  we  see  that  the 
AF  and  the  WD  are  related  by  a  transform  sim.ilar  to  a 
2-dimensional  F.T.  (Figure  14). 


yrr(u,v)^^ 

X(u,t)-  W(t,v) 

^  •Y<t,T)^ 


Figure  14 


Thus 

(1) 

W(t,v) 

_1 

IJ  X(u,T)e^^^^  '^^^dudx. 

and 

(2) 

X  (u  ,T  ) 

_1 

1  'W(t,v)e  ^  ^  ^dtdv 

There  are  three  basic  differences  between  the  WT 
and  the  AF.  First,  the  WD  is  always  real-valued,  while 
the  AF  is,  in  general,  a  complex-valued  function.  The 


modulus  is  frequently  used  in  applications  so  it  is 
reasonable,  in  many  cases,  to  compare  the  miodulus  of  the 
AF  to  the  VJD.  Only  when  f  is  real- valued  and 
symmetric,  is  the  AF  the  samie  as  the  WD  (up  to  a 
scale  factor) .  This  was  demonstrated  in  Examples  1  and  10 

The  second  difference  between  the  AF  and  the  WD 
is  the  manner  in  which  these  functions  transform  trans¬ 
lated  functions.  The  modulus  of  the  AF  is  invariant  to 
translations  of  either  f  or  its  F.T.  The  WD  shifts 
as  its  generating  function  shifts.  For  this  reason, 
Claasen  and  Mecklenbrauker  conclude  the  AF  is  not  well 
suited  for  general  signal  analysis  [13]. 

Finally,  the  AF  and  VJD  each  characterize  signal 
duration  and  bandwidth  of  the  generating  function  f  in 
different  ways.  Effective  duration  and  bandwidth  of  f 
are  equivalent  to  partial  derivatives  of  the  AF  but 
equivalent  to  the  global  moments  of  the  WD.  Consequently 
the  WD  was  shown  to  characterize  directly  the  Heisenberg 
uncertainty  principle.  On  the  other  hand,  the  AF  has 
an  associated  uncertainty  principle  (property  11)  which 
is  an  analogy  to  the  Heisenberg  uncertainty  principle. 

The  WD  and  AF  have  three  basic  similarities. 
First,  each  time- frequency  transform  of  f  is  closely 

A. 

related  to  the  t-f  transform  of  F  =  f.  The  t-f 
transform  of  F  is  a  simple  rotation  of  the  t-f  trans¬ 
form  of  f.  As  such,  Hermite  polynomials  are  sometim.es 
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associated  with  these  functions.  For  instance,  an 

2 

orthogonal  system  on  L  based  on  Hermite  polynomials 
was  used  by  Klauder  to  achieve  highly  concentrated  AFs 
about  the  origin  [8],  This  same  orthogonal  system  was 
used  by  De  Bruijn  in  the  proofs  of  his  inequalities  for 
the  V.’D  C  5  ]  . 

Another  similarity  is  that  both  functions  Transform 
quadratic  operations  on  the  generating  function  to  linear 
operations  on  the  AF  and  WD.  Therefore,  the  decompo¬ 
sitions  of  the  AF,  introduced  in  Section  7.2,  may  be 
applied  to  the  VJD. 

Finally,  both  t-f  transforms  yield  similar  constants 

when  integrated  over  the  entire  plane.  In  this  case  we 
2 

compare  |x|  to  the  WD.  For  clarity  we  repeat  (4.19) 
and  (9.16). 

(4.19)  —•  |||X(u,t)  I  ^  dudT  =  llfll^, 

and 

(9.16)  jj  W(t,v)dtdv  =  'Ifl:^. 

2 

Let  us  define  a  norm  on  IR  as 

(3)  I'Xl^  =  (2“  |||X(u,T)|  ^  dudT)^^^. 

2  2  2 

Then  X  ^  L  ( P,  )  because  f  €  L  (P).  Moyal's  formula 
is  [3], 


(4) 


l!f!!  . 


—  I  W  (t ,v)dtdv  = 

2tt  j  j 

2 

This  gives  us  a  Parseva] -like  relationship  in  IR  ; 


2  2 

Therefore,  in  the  space  of  L  ( IR  )  ,  with  the  defined 
norn,  the  transform  (13.1)  is  an  isometry. 

§14.  Conclusions 

Time- frequency  functions  like  the  AF  and  WD  are 
imiportant  in  pure  mathematics  and  engineering.  They  are 
functions  which  transform  a  function  of  time  into  a  func¬ 
tion  of  time  and  frequency.  They  inherently  embody  much 
of  the  theory  of  Fourier  transforms.  Study  of  these 
functions  provides  a  richer  understanding  of  the  uncertainty 
principle  and  is,  therefore,  valuable  in  the  field  of 
harmonic  analysis.  In  the  applied  fields,  a  time-frequency 
transform  helps  to  visualize  the  frequency  content  of  non¬ 
stationary  functions.  These  transforms  help  explain  the 
intimate  relationship  between  the  time  and  frequency 
components  of  a  signal. 

The  AF  is  a  tool  routinely  used  in  radar  waveform 
design.  It  is  also  used  in  developing  the  theory  of 
Fresnel  diffraction  and  Fourier  optics  [13].  It  was  shown 
to  provide  a  better  understanding  of  linear  frequency 
modulation  and  pulse  compression.  The  AF  is  an  important 
function  in  both  the  applied  and  theoretic  fields. 


The  \'!D  has  been  for  many  years  an  asset  to  the 
field  of  quantum  physics.  Recent  work  has  shown  it  to 
be  a  useful  tool  in  signal  analysis  and  suitable  for 
hardware  implementation  [3].  It  was  also  shown  to 
generalize  two  common  notions  in  signal  analysis; 
instantaneous  frequency  and  group  delay.  The  WD  is 
incompressible  and  reflects  the  Heisenberg  uncertainty 
principle .  The  WD  is  also  an  important  function  in 
theoretic  and  applied  fields. 

Time-frequency  transforms,  like  those  presented  in 
this  paper,  have  a  broad  application.  It  is  the  opinion 
of  this  author  that  such  functions  may  have  even  more  to 
contribute  in  such  fields  as  spectrum,  estimation,  and  are 
worthy  of  future  study. 
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Appendix  A .  The  Complex  Signal 

Frequent  reference  is  made  throuphout  this  paper  to 
the  complex  representation  of  a  real  signal.  This  is  also 
known  as  the  Gabor  representation  and  the  analytic  signal. 
In  this  appendix  we  will  define  the  analytic  signal  and 
explain  v/hy  it  is  a  useful  tool. 

In  most  applications  of  signal  processing  techniques, 
the  signal  of  interest  is  a  real-valued  function  of  time, 
say  s(t).  As  was  stated  in  (3,1),  for  radar,  s  may  be 
modeled  as 

(1)  s(t)  =  a(t)cos(a)Qt+(})  (t ) )  . 

This  description  is  mathematically  cumbersome  and  its 
Fourier  transform  has  a  lot  of  redundant  information  since 

(2)  S(-w)  =  STO. 

ConsequeiiLly  a  complex-valued  function  f  is  formed  so 
that  Re{f}  =  s.  The  Gabor  representation  is  such  a 
com.plex- valued  function  and  is  defined 

f(t)  =  s(t)  +  is(t) 

V  V 

where  s  is  the  Hilbert  transform;  of  s.  s  is  defined 

(4)  s(t)  =  -  )  — ~  dr  =  s(t)  * 

TT  r-t  mt 

The  Hilbert  transform  has  the  useful  property  that 
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(5)  S(co)  =  -i  sgn(w)S(w)  , 

V  V 

where  S  is  the  Fourier  transform  of  s. 

This  mav  be  verified  by  lettine  x(t)  =  so  that 

TTt 

s(t)  =  s(t)  *x(t) 

and 

SCw)  =  S(to)X(u)). 


Then  (A.  5)  follovjs  from  the  Fourier  transform  pair 
1 

^ — ►  “I  s^n  u'j  • 

TTt 

From  (A.  3)  we  see  that 


{2S (to )  to  >  0 

3(0)  to  =  0 

0  to  <  0  . 


The  Gabor  representation  solves  the  problem  of  redundancy, 
but  using  (A. 3)  in  calculations  may  be  equally  as 
cumbersome  as  using  the  original  s(t).  Consequently,  in 
many  cases  an  approximation  to  (A. 3)  is  used  in  calculations: 

i  (to  t  +  (})  ( t )  ) 

(7)  f (t)  =  a(t )e 

This  approximation  becomes  exact  for  narrowband  signals. 

That  is,  define 

F^(u))  =  I  I  f(t)e"^'‘''^  dt 

and 

F_(u))  =  i  j  dt. 


Then  f  is  narrowband  if 


F_|_(u3)  =  0  for  all  u)  <  0 , 

and 

FCuj)  =  0  for  all  co>0. 


In  general,  the  error  in  the  approximation  (A. 7)  is 


e (t)  =  21 


m 


r  1  rO 

( 

J  -oo 


(  ..  iu)t  , 
(w;e  dcoy. 


This  is  obviously  zero  for  narrowband  signals. 
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ArT:^endix  B.  The  Doppler  Appro xijr.ati_on 

This  appendix  will  explain  the  doppler  approximation 
and  why  the  effect  of  a  moving  target  on  a  narrovjband 
signal  may  he  modeled  as  only  a  shift  in  the  carrier 
frequency . 

A  radar  transrr.its  an  electrom.a met i r  sirnal,  vihich 
v;hen  striking  suitable  surfaces  is  reflected  and  refracted 
sim.ilar  to  light.  If  the  surface  is  moving,  the  reflected 
frequency  of  the  signal  v;ill  appear  to  be  different  than 
that  which  v;as  transmitted.  This  is  analogous  to  the 
comiiaonly  deserved  phenomena  of  a  constant-pitch  train 
whistle  appearing  high  as  the  train  approaches  and  low  as 
the  train  passes.  This  phenomena  is  known  as  the  doppler 


effect . 

In  the  case  of  radar,  should  the  target  be  stationar;, 
(in  the  sense  that  the  target  velocity  vector  has  no 
component  in  a  radial  direction  to  the  radar) ,  then  the 
tim.e  delay  of  the  returned  signal  will  be 


where  P  is  the  constant  range  and  c  is  the  velocity 
of  propagation.  The  returned  signal  is  then  a  time  delayed, 
suitably  attenuated  version  of  the  transmitted  signal.  If 
s  (t)  is  the  transmitted  signal  then  the  returned  signal  is 


S  (  t  )  =  S_^  (t-T  )  . 

r  t 


If  the  reflecting  surface  is  moving,  then  the  rang 
becomes  a  function  of  time  and  therefore  the  time  delay 
becomes  a  function  of  time.  Then 

s  ( t )  =  s ,  ( t-T ( t ) )  . 

r  t 

The  exact  relation  between  the  time  delay  T(t)  and  the 
range  to  the  target  is  [6] 

(1)  T(t)  =  -  R[t  -  T(t)  ]. 

c  2 

Defining  T(t)  =  t,  (B.l)  can  be  expanded  in  a  Taylor 
series  about  this  point; 

T(t)  =  T(t)  +  f(T)(t-T)  +  I  f(T)(t-t)^  +  ...  . 

The  returned  signal  is  then 

S^(t)  =  S^ ( t-T( T )-f( T ) ( t- T )  f(T)(t-T)^  -  ...). 

Use  T(t)  =  T  and  this  is 

(2)  S  (t)  =  s([l-t(T)][t-T]-'T(T)[t-T]^  -  ...). 

r  t  2 

Use  Standard  notation  for  velocity  and  acceleration,  i. 
v(t)  =  R(t)  and  a(t)  =  R(t),  and  differentiate  (E.l) 


to  see  that 


T(t  ) 


(U) 


[l+(-^)v(l)  ]■ 

c  2 


In  radar,  the  velocity  of  propaFation.  c,  is  nearly 
the  speed  of  ]ight,  therefore  c  >>  v(t).  (B.5)  and  (5.4 

are 

(5)  T(t)  =  (-)v(^) 

c  2 

(6)  T(t)  =  (-)a(I). 

c  2 

By  definition,  the  doppler  effect  is  the  linear 
stretching  of  the  time  variable  in  (B.2).  Therefore  we 
only  consider  the  first  term  in  the  series  (B.2).  This 
is  equivalent  to  assum.ing  the  target  velocity  is  constant 
near  t  so  that  v(t)  =  v(^-)  .  The  linear  approximation 
of  (B.2)  is 


(7)  s  (t)  =  s^[w(t-T)], 

r  t 

where 

.  2 

w  =  1-T(t)  =  1--v(t). 

c 


Finally,  assume  s^ 
we  may  write,  for  fixed 


is  a  narrowband  signal. 

w , 


Then 


(3) 


icot 

s  (t)  =  u(t)e 

t 


v;here 
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Then  the  doppler  stretching  vj,  of  the  time  variable 
of  the  returned  signal  may  be  assumed  to  apply  only  to 
the  carrier,  This  seems  justified  for  the 

perturbation  of  t  vjould  be  less  perceptible  on  the 
slowly  varying  function  a(t)  and  G(t)  relative  tc  the 


fast  varying  carrier.  Hence 


s  ( t ) 
r 


u( t-T  )e 


icj(w)  (t-T  ) 


Define  d  =  bo  be  the  doppler  shift  then 


s  (t)  =  u(t-T)e 

r 


i(cj-4>)(t-T) 
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